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ABSTRACT 

Collisionless relativistic shocks have been the focus of intense theoretical and numerical in- 
vestigations and these interesting physics have a direct impact on the generation of energetic 
particles and the interpretation of gamma ray spectra. The Fermi acceleration process that 
takes place in these shocks is intimately linked with the excitation of micro-turbulence respon- 
sible for the shock formation, electron heating and supra-thermal tail generation that in turn 
excites micro-turbulence, developing thus a self-sustaining phenomenon. In this paper we dis- 
cuss the development of the micro-turbulence and we investigate two important issues: firstly 
the transport of supra-thermal particles in the excited microturbulence upstream of the shock 
and its consequences for the efficiency of the Fermi process; secondly, the preheating process 
of the incoming background electrons as they cross the shock precursor and experience rela- 
tivistic oscillations in the electric field of the micro-turbulence. We emphasize the importance 
of the motion of the electromagnetic disturbances relatively to the background plasma and 
to the shock front. The investigation is carried out for the two major instabilities involved in 
the precursor of relativistic shocks, the filamentation instability and the oblique two stream 
instability. Bearing in mind these new results, we analyze the performance of the Fermi ac- 
celeration process in various high energy sources and especially in the termination shock of 
gamma-ray bursts. We emphasize the high efficiency of a collisionless ultra-relativistic shock 
for accelerating electrons and radiating a synchrotron-like photon spectrum up to several GeV. 
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1 INTRODUCTION 

The micro-physics of collisionless relativistic shocks has been in- 
tensively investigated during these last years, both through numer- 
ical wor ks and theo r etical o nes, especially since the remarkable 
work of ISpitkovskvl j2008al lbh. The physics of those shocks de- 
scribes the interplay between three phenomena, namely, the shock 
formation by reflecting back a fraction of the incoming particles, 
the generation of micro-turbulence by the reflected or back scat- 
tered particles that interact with the incoming plasma, and the de- 
velopment of the Fermi process. Much progress has been achieved 
within a few years through the development of theoretical and nu- 
merical works. Although the numerical simulations cannot reach 
the highest energies involved in the astrophysical phenomena, these 
works are already precise enough for inferring astrophysical con- 
sequences. 

The physics of a relativistic shock is governed by few param- 
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eters, its Lorentz factor T s , the density of the ambient or upstream 
medium n and its magnetic field B . These two parameters deter- 
mine two essential physical quantities, namely, the mass energy 
density pc 2 which fixes the incoming energy to the shock front, 
T 2 pc 2 , and the magnetization parameter a that measures the ratio 
of magnetic energy flux crossing the shock over the incoming en- 
ergy density: 



B 



AnT]pc 2 Anpc 2 



sin 2 6»b , 



(1) 



where fi,| f = T S B sin 6 B is the transverse component of the back- 
ground magnetic field in the shock front rest frame, B being the 
angle of the magnetic field with respect to shock normal, measured 
in the upstream rest-frame. Numerical simulations provide two es- 
sential parameters for astrophysical applications, namely the con- 
version factor £ cr of the incoming energy into cosmic rays, and the 
conversion factor £b into magnetic energy: 



P a = ^Tlpc 2 
An = ^ PC 



(2) 
(3) 
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where the cosmic ray pressure P cr and the level of magnetic turbu- 
lence S|f are measured at the shock front. These two crucial param- 
eters £ CI and£s are expected to be on the order of 1 - 10 percent 
dSironi & Spitkovsky||201l3) . Actually the cosmic rays are consid- 
ered, as shown by numerical simulations and explained by theory, 
as the source of magnetic, and more generally electro-magnetic, 
turbulence. Indeed the collisionless shock forms by reflecting back 
a fraction of incoming particles on a barrier that can be due to the 
mean field, helped by an electrostatic potential when the plasma 
in composed of protons and electrons, or by the growth of elec- 
tromagnetic waves, that in turn, triggers micro-instabilities. Parti- 
cle acceleration has been observed in PIC simulations of unmag- 
netized relativistic shocks ( Spitkovsky 2008b; Keshe t et al. I l2009l ; 
iMartins et ai]|2009l ; ISironi & Spitkovsky||2009i. l2011ah . and indeed 
one must expect the develop ment of the Fermi process when the 
magnetization is very weak dLemoine & Pelletier! l20~ioh , because 
microturbulence can then grow and provide the necessary scatter- 
in g; this picture has been valid ated by the recent PIC simulations 
of Sironi & Spitkovskvl ( 12011 j) . In this context of a proton elec- 
tron plasma of low magnetization, the reference time scale is 0)~l 
and the spatial scale of reference is the inertial scale of protons 
8, = c/co pi . 

A collisionless relativistic shock at low magnetization is 
mostly governed by the development of two micro-instabilities 
that take place in the precursor, generated by the reflection of a 
part of inc oming particles, namely the Weibel or filamentation 
instability dMedvedev & Loeblll999h and the oblique two stream 
instability (hereafter OTSI; lFainberg et alj[l970l. iBret et alj|2005n . 
Because the upstream mean magnetic field in the front frame is 
generically almost perpendicular to the shock normal, it limits 
the penetration of reflected particles in the upstream flow; even 
weak it puts severe constraints to the growth of the instabilities 
by limiting the interaction length and the higher the Lorentz 
factor of the shock the weaker the value of the magnetization 
paramete r below which the growth of th e unstable modes is 
allowed dLemoine & Pelletier|[2oTol 1201 lah . In the presence of a 
mean field, the efolding length scale (f g ~ £^i l2 S\) is shorter than 
the length of the foot, or the penetration length of reflected protons, 
tfoot = LoU -/3 s )/r s ~ 6J(2T S s/cr), when the magnetization is 
as small as a < ^ a /4I^. The most unstable mode with decreasing 
magnetization is th e OTSI with the condition <x < f c / 3 /p 1 ' 3 !^, 
where n = m e /m p dLemoine & Pelletie3l2O10h . Its role in the de- 
velopment of relativistic shocks is still unclear and this paper aims 
at shedding light on this iss ue. Then for weaker magn etization, 
precisely when a < fcr/Tj dLemoine & Pelletierll2010h . Weibel 
modes can grow and the turbulent state so created allow particle 
scattering off small scale magnetic disturbances (filaments), and 
relativistic Fermi acceleration can develop. The scattering and 
diffusion of the particles in a small scale static and isotropic 
turbulence superposed on a mean field has be en investigated 
nume rically and theoretically in a previous paper dPlotnikov et al .1 
1201 lb . In this paper we emphasize the issue of the motions of the 
magnetic disturbances and their anisotropy. With the studies of 
particle transport in both types of micro-turbulence, Weibel and 
OTSI, we draw their consequences for the collisionless relativistic 
shock formation, for electron pre-heating and for the performance 
of the Fermi process at astrophysical relativistic shocks. 

For completeness, we mention that there are also relevant in- 
stabilities triggered by the compensation current in the background 
plasma that compensates the current carried by the reflected par- 
ticles either along the mean field (but not generic) or across the 



mean field (unavoidable). For instance, reflected particles of oppo- 
site charge rotate in opposite direction in the transverse mean field 
and thus produce a very intense diamagnetic current responsible 
for a Buneman instability dLemoine & Pelletierl2011ah . and the in- 
stability is rapidly growing. Those current instabilities produce a 
turbulent heating of the electrons up to some temperature that re- 
duces the anisotropy of the electron distribution function so that the 
instability saturates. Therefore these current instabilities participate 
in the process of preheating electrons, that arrive at the shock front 
with a relativistic temperature. 

A collisionless relativistic shock can be described as a self- 
sustaining structure in which reflected particles generate a grow- 
ing turbulence level that is responsible for particle reflection, for 
particle scattering and particle heating. Thus the key point is the 
identification of the relevant instabilities (one or two main) that 
produce modes able to accomplish this task. The main focus is on 
the Weibel instability that creates magnetic fluctuations on micro- 
scales, which stand in a frame that has to be specified and the OTSI 
that creates electromagnetic waves with an important electrostatic 
component in the background plasma. The waves excited in the 
background plasma cannot carry an energy density exceeding its 
mass energy density, so it can be written t/ w = <f w pc 2 with £ w < 1. 
In the frame moving with the shock front, the energy density is 
r^ w pc 2 , the conversion factor into waves f w being expected to be 
a sizable fraction in order to participate at both the shock formation 
and the generation of a supra-thermal tail. The pressure of relativis- 
tic particles can reach a comparable level, even a little larger since 
they are the cause of the instabilities. The first generation of re- 
flected particles constitutes the main content of supra-thermal par- 
ticles that pervade the ambient plasma with an energy density much 
larger than pc 1 as measured in the frame of the ambient plasma. 
Its interaction with the background plasma can be described per- 
turbatively as long as the proton beam of (apparent) density «b 
and Lorentz factor yb is such that <jj p \, <k a> p , which amounts to 
iib/jb *K n/ji; now for a beam reflected by a shock, n b = ^^T^n and 
7b ~ Tl; thus the weak interaction criterium becomes £ cl <K 1//j, 
which is always realized. 

The weak interaction of the very energetic beam with the 
ambient (or upstream) plasma renders the calculation of instabil- 
ity convenient in the frame of the background plasma, but, be- 
cause of the tremendous energy of the beam, a saturation process 
through beam particle trapping or resonance broadening cannot be 
expected, as will be checked further on. The only nonlinear effects 
that can be relevant in the context of relativistic shock physics are 
nonlinear wave interactions, in particular the effects caused by the 
ponderomotive force associated with the instability wave. 

The reflection of incoming protons, electron preheating and 
the motion of scattering centers are entangled issues and will be in- 
vestigated step by step. The paper is organized as follows. In Sec- 
tion 2, we emphasize the important issue of the motion of the frame 
in which the Weibel filaments are static, an issue which has not 
been addressed yet, to our knowledge. We also summarize previous 
findings on a similar issue for the OTSI mode. We investigate the 
influence of the motion of the electromagnetic modes on the reflec- 
tion process at the shock front. In Section 3 we study the transport 
of supra-thermal particles in both Weibel and OTSI turbulence, on 
the basis of numerical simulations of test particle propagation. Sec- 
tion 4 is devoted to the issue of electron pre-heating. We show that 
the relativistic oscillation of the incoming background electrons in 
the electric field of the micro-turbulence modes lead to efficient 
preheating on a short timescale. In Section 5 we apply our results 
to the concrete case of the external relativistic shock of a gamma- 
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ray burst; we emphasize their capability to efficiently accelerate 
electrons so that it can radiate an extended synchrotron-like spec- 
trum up to GeV photons. However we briefly show the difficulty 
to achieve ultra-high energies for protons in those shocks. In Sec- 
tion 6 we discuss possible sites of ultra-high energy cosmic rays 
(UHECR) generation. We summarize our results in Section 7. 



2 MICRO-TURBULENCE IN THE PRECURSOR OF 
COLLISIONLESS WEAKLY MAGNETIZED 
RELATIVISTIC SHOCKS 

2.1 The frame of magnetic filaments 

Consider first the growth of magnetic perturbations due to Weibel 
instability triggered by a beam of velocity u b and Lorentz factor 
7b = (1 - v^/c 1 ) ~ 1/2 ~ r 2 (upstream frame), interacting with a 
background plasma of density n at rest; at a shock wave, the beam 
of returning particles carries an energy density (; a y^nm p c 2 , so that 
the beam plasma frequency a) pb ~ ^ r '/ 2 / j 1/2 ai De , with qi pe the back- 
ground electron plasma frequency (see lLemoine & PelletierJ [20 1 dl 
for details). In the upstream rest frame, the efolding length scale of 
the instability is written 



the detailed growth rate (see appendix A) being 

Pbk±S e 



7mst = 3^ = Vicr" 



(1 + k 2 J 2 ) 1 ' 2 



(4) 



(5) 



In the above expression and throughout this paper, k ± represents 
the wavenumber component transverse to the shock normal, i.e. 
tangential to the shock front, while the (longitudinal) wavenum- 
ber component along the shock normal is written k^. For the fila- 
mentation mode, k ± » k\\. In general, one takes the limit k\\ — > 0, 
which leads to the above aperiodic mode with %a> = 0. How- 
ever, for a small, but finite longitudinal wave number k\\, the Weibel 
modes have a non-vanishing real frequency which indicates that 
these magnetic filaments propagate at high velocity. Considering 
first the case of a cold electron background: 



(O r = k\\Vb 1 -gcrM 



1 + k 2 6l 



(6) 



We remark that this phase velocity is consistent with the result ob- 
tained in the center of mass frame of the counter streaming config- 
uration when 7b - 1; for /j = 1 and f CT = 1/2, we find = u b /2, 
which corresponds to a vanishing phase velocity in the center of 
mass frame. In the upstream rest frame, these waves are quasi- 
electromagnetic waves with the electric field E mostly oriented 
along the beam (shock normal), a dominant magnetic component 
perpendicular to both k ± and E, and finally a small electrostatic 
component of E along k ± . The magnetic component remains the 
most intense component despite a phase velocity close to c along 
the normal direction because B = (k ± /k\\)E, and the energy density 
is mostly magnetic. The fact that they have a phase velocity close 
to c implies that they suffer negligible Landau damping (as long as 
the thermal velocity of the electrons is small). 

If the background electrons have been preheated to relativis- 
tic temperatures, their electromagnetic response and therefore the 
instability are modified. The maximum growth rate remains un- 
changed, 3<y ~ £a 2 (i> p i, but the spatial scale at which maximum 
growth occurs is now larger because the electron inertial scale is en- 
larged by a factor with y e the mean electron thermal Lorentz 



factor in the upstream rest frame y e = 1 + 3T e lm e c 2 . This lat- 
ter spatial scale tends to 6\ in the limit of equipartition, meaning 
y e — > m p /m„. In a relativistically hot background, the phase veloc- 
ity of the modes along the beam direction is considerably slowed 
down, but remains nonetheless relativistic for reasonable values. If 
one defines y m as the Lorentz factor associated to the mode velocity 
along the beam direction, the calculations presented in Appendix A 
indicate 

1-1/2 



7m - \(%rf a ) 1 ' 2 + 2fo 



(7) 



Equipartition corresponds to y e fi = 1, and prior to equipartition 
y e /i < 1. For relativistically hot electrons far from equipartion, 
7m - (2&r)~ 1/2 , then y m ~ (^fife)' 1 ' 4 when y er i > and finally 
the lowest value of y m ~ 2 is reached at equipartion. This value is 
much smaller than that obtained in the cold electron background, 
7m ~ (^civ") l/2 , but it nevertheless plays an important role in what 
follows. 

Note that the filamentation instability in the usptream of a rel- 
ativistic shock is quenched when T s < (f cl /u)~ 1/2 as long as the 
electro ns are kept cold due to the finite angular dispers ion of the 
beam dRabinak et~al]|201 ll iLemoine & Pelletiej|201 lah . In a hot 
background however, the filamentation instability remains strong, 
even when the beam angular dispersion is taken into account. 

So clearly, the Weibel filaments that can be considered as wave 
packets of transverse size £ ± ~ y l J 2 S e and extension in the normal 
direction fy > y l J 2 S e , but necessarily finite (at least because of the 
finite growth length in a precursor of finite extent), are unavoidably 
moving at high speed with respect to the ambient medium. This 
is a very important point that has not been taken into account in 
the literature, to our knowledge. They move at a highly relativistic 
velocity when the electrons are cold, with a Lorentz factor y m = 
(£ crt u)~ 1/2 , which can be of several hundreds (typically for £ cr = 
0.1, 7m - 140). When the electrons are relativistically hot, y m is 
more moderate. Three cases deserve attention: y m ~ T s where the 
filaments move more or less at the same speed as the shock front, 
7 m < r s where the shock front catches up the filaments and the case 
where the filaments can even run faster than the shock for T s < y m . 
In this latter case the generation of filaments might lead to shock 
reformation. 

The filaments are static in the frame moving at B m with re- 
spect to the background plasma, however in this frame they are no 
longer purely magnetic and comprise an electrostatic field of al- 
most the same intensity as the magnetic field: E' ± ^ y m B = B', 
since B' = y m B, with the prime denoting the quantity in the 
filament rest frame. The electromagnetic energy density in this 
frame is thus U' em = y 2 n B 2 /(4n). The electromagnetic flux that 
crosses the shock front, measured in the co-moving front frame is 

5e.m.tf = V s r 2 fi 2 /47T. 

The relative velocity of the filaments with respect of the shock 
front can be written 



An|f 



r, 



2 -r 2 



7 2 +r 2 

/ m 1 A s 



and the associated Lorentz factor 



2 \ T s y m 



(8) 



(9) 



Let us look at their growth rate as measured in the shock front 
frame. As shown in the appendix B, the product 7 e 7i nst is invariant 
under Lorentz transformations, where 7 5 is the Lorentz factor as- 
sociated with the group velocity. Thus the transformed growth rate 
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7insl|f 



V2a 



r s y m 

' T""9 , 9 ' 
r s + 7m 



(10) 



which is maximum for y ra = r s , when the magnetic disturbances 
are standing in front frame. 

The motion of the wave packets carries special importance 
when determining the condition for particle reflection, the level of 
electron pre-heating and also for the analysis of the Fermi process, 
as discussed in the following. 

In the case of an ultra-relativistic shock in an electron positron 
plasma, y m is at most a few and the shock front catches up the 
magnetic disturbances (T 8 » y m ). The behavior is the same in a 
proton-electron plasma when 1 <k y e < m p lm e . 

The discussion about the motion of magnetic disturbances 
is also important for their transmission in the downstream flow. 
Indeed the relative motion of the upstream flow with respect to 
the downstream one is characterized by a relative Lorentz factor 
T r = r s / V2. Because the modes propagates at a speed u m with 
respect of the upstream plasma, the Lorentz transform of their fre- 
quency and normal wave number to the downstream frame leads 
to 

co' = r r (/3 m - B Y )k z c - ^T T k z c ( p - -7 



1 / 1 1 

k' z = r r (i-Aj8 m )i^ + - 



(11) 



Thus, when y m <k T s , the modes are perceived as electromagnetic 
waves of the vacuum in the downstream frame, propagating back- 
wards (a>' -k'.c). In the opposite extreme, they are perceived as 
forward propagating vacuum waves. The interesting case is when 
Jm ~ T s , where they are perceived as low frequency magnetic per- 
turbations. However, one should stress here that the motion of the 
filament is calculated at the linear level, while the transmission of 
modes downstream proceeds in the highly nonlinear regime. There- 
fore it is not clear at present whether the above Eg. II llapplies to the 
transition; it should be taken with caution. 



2.2 The role of OTSI turbulence 

A priori the Weibel instability can produce all the effects expected 
to occur at a relativistic shock including Fermi acceleration. How- 
ever the OTSI is unavoidable, and it even grows a little faster than 
the Weibel instability. The OTSI is a resonant instability of elec- 
tron plasma waves (a> ^ ai p( ,) with a sharp selection of the wave- 
vector component along the beam: k- = co pe /c, and that grows much 
faster than the usual two stream instability (which has no trans- 
verse component of its wave-vector) when the tran sverse compo- 
nent of the wave-vecto r is of order co pe /c also (see iFainberg et all 
1 197(1 lBretetal]|2005l) . These modes are mostly electrostatic in 
the background plasma frame and partly electromagnetic and grow 
a little fas ter than Weibel mode in the precursor of a relativis- 
tic shock dLemoine & Pelletie3l2010h , their frequency is slightly 
shifted by the beam with to = a> pi ,(l - \S\/2) and their growth rate 
3co = V3 2- 4/3 |<5|w pf with \6\ = fc/V /3 > whereas the growth rate 



of the Weibel modes scales as f cr ' n 



l l 2 ,M 2 , 



in the same cold back- 



ground limit. The ratio of the electromagnetic component over the 

electrostatic one is of order |<5|. 

As discus s ed in iLemoine & Pelletieil d2011al) . 
IShaisultanov et al. the two stream instability becomes 

inhibited once the electrons are heated to ultra-relativistic tem- 
peratures and the Weibel instability becomes the dominant 
mode. 



In order to study the dynamics of particles in the modes, it 
is interesting to move to the wave-frame where the particles ex- 
perience a static electric field and a static magnetic field. The 
wave frame velocity with respect to the background plasma is 
B m = 1 - \S\I2 and the corresponding Lorentz factor 

i<r 1/2 : 



,.-1/6 1/6 



(12) 



with respect of the background plasma. These mildly relativistic 
wave-packets are thus rapidly overtaken by the shock front. 

As in the case of Weibel modes, the dynamics in OTSI modes 
is also governed by a couple of transverse fields of similar ampli- 
tude when one shifts to the rest frame in which these modes are 
static. A crucial difference however is that the OTSI mode appears 
as a high frequency quasi vacuum wave in the front frame, with a 
unique wavenumber along the normal direction, revealing a clear 
periodic pattern in that direction. 



2.3 The reflection of incoming particles 

The deviation of the incoming protons becomes significant when 
the level of magnetic fluctuations is such that the effective Larmor 
radius Fl ~ where 6\\ denotes the coherence length of the mag- 
netic disturbance in the normal direction. In order to properly state 
that condition, we must express it in the co-moving frame of the 
magnetic disturbances; indeed if one moves to another frame ac- 
cording the Lorentz transformation, the disturbances become time 
dependent and fully electromagnetic. Therefore we state that re- 
flection is achieved when, in the proper frame of magnetic distur- 
bances, r L | ra ~ ^|| ,„; that estimate amounts to eB\ m £\y m ~ y m m p c 2 
with = y m {\\ and B\ m = y m B. Since the electric field is of sim- 
ilar amplitude in the filament frame, the corresponding condition 
with the electric field component is that the energy of the incoming 
protons is comparable with the work done by the electric force over 
the scale €«. Therefore 



An 



5\ 



(13) 



which leads to a magnetic conversion factor, as defined in the in- 
troduction, 



1 6, 



Y m <"|i <"|i 



(14) 



the second estimate stems from the fact that the reflection occurs 
in the phase where the electrons are relativistically hot and thus 
7m ~ (fcrTe/")" 1 ' 4 an d the last one when y e fi < £ cr . Furthermore we 
expect that the wave extension £\\ (defined in the usptream frame) is 
not larger than the growth scale i g = ^.^ 2 S\, so that mode envelopes 
on scale of a few 6, can produce the necessary reflection. 

Under that condition a significant amount of incoming pro- 
tons are reflected back and self-sustain that magnetic turbulence 
through the Weibel instability. When those conditions are fulfilled, 
the self-sustaining turbulence through the reflection of a fraction of 
incoming protons stabilizes at a level that remains much smaller 
than all other processes of nonlinear saturation such as proton trap- 
ping, resonance broadening. Indeed, the trapping of reflected beam 
protons requires that the level of micro-turbulence is such that 
the time scale of non-linear oscillation r m = 4(£/£*)' /2 / c (where 



eE'{ c ) reaches the growth time w , ; since w p i 



/.ill 



this amounts to ~ f CT 6^/(5?, which is totally impossible with 
beam protons of energy e ~ T 2 m p c 2 . A similar conclusion can be 
drawn when one considers the resonance broadening effect, which 
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stems from the fact that particle scattering or diffusion broadens 
the resonance of modes with the beam responsible for the instabil- 
ity. Indeed v s > toi p b would require > £c/ 4 (/ c /(5i) 1/2 f, which is 
impossible with a so high energy of beam protons. 



2.4 Heating and transport in the micro-turbulence 

In the Section that follows, we discuss in some details the scattering 
of supra-thermal particles and the heating process of background 
particles in the motional micro-turbulence upstream of the shock. 
These actually represent two different facets of a similar problem, 
namely particle transport in a time varying micro-turbulence. How- 
ever, in the test particle picture that we adopt in the following, these 
two populations, the supra-thermal particles and the background 
electrons, differ one from the other by their wiggler parameter 



eE'l\ 



(15) 



as expressed in terms of the micro-turbulent electric field E' and 
transverse scale C' ± in the wave frame. As discussed in the follow- 
ing, in the wave frame E' and B' are of the same order, so that there 
is no ambiguity in the definition of a. 

Cold background electrons have a Lorentz factor y' ~ y m in 
the wave frame, so that 



y' m e d\ 



(16) 



meaning that these electrons experience relativistic oscillations on 
a coherence length scale in the wave frame. Of course, as the elec- 
trons near equipartition, the ratio a/y' becomes closer to unity. 

Consider now a supra-thermal electron, with Lorentz factor y 
in the upstream rest frame, becoming y' = y/y m in the mode rest 
frame (assuming y » y m ). The minimal Lorentz factor of supra- 
thermal electrons is y min ~ T 2 m p /m e in the upstream frame, hence 
one can write 



2_ „ yl/2 7m 7min l± 

Y ~ Tj y 6\ ' 



(17) 



which is expected to be much smaller than unity: recall that £ x is 
expected of the order of 6 e if the background electrons are cold (in 
which case y m can be large if the Weibel instability is not quenched 
by the angular dispersion of the beam), but of order i5j if the back- 
ground electrons reach equipartition with the ions, in which case y m 
becomes of order of a few for the Weibel modes. Thus a/y' <k 1 
for supra-thermal electrons, while a/y' » 1 for the background 
electrons, typically. 

Before discussing the transport of these particles in the up- 
stream micro-turbulence, it may be useful to briefly recall how 
transport is thought to take place downstream of the shock. Accord- 
ing to the numerical simulations of Chang et al. (2008), the micro- 
turbulence downstream of the shock ap pears mostly stat i c and mag- 
netic in the downstream rest frame. In lPlotnikov et alj ( l201ll) . we 
have discussed in some detail the transport coefficients of high en- 
ergy particles in such a micro-turbulence, accounting for the pres- 
ence of a transverse background magnetic field. 

One first defines the scattering frequency 



2 c eB£ c 
34 



(18) 



with £ c the coherence scale of the fluctuations, B = ■JB^ + SB 1 
the rms magnetic field and e the particle energy, in the downstream 



frame. When there is no mean field, the spatial diffusion is isotropic 
and is simply described by the following diffusion coefficient 



D = 



I EI 

3 v s 



\eB£ c 



(19) 



When there is a mean field, even weak, the spatial diffusion is 
anisotropic, because the scattering frequency becomes smaller than 
the Larmor frequency in the mean field at high energy, meaning that 
e » eBC c . The diffusion along the mean field, D\\ remains similar 
to the previous one, but the transverse diffusion D ± is given by 



(20) 



with o) L0 = eB C c /€, similarly to classical diffusion. Thus instead 
of increasing like e 2 at high energies, it reaches a maximum finite 
value 2c£ c B 2 /(9B 2 Q ). 

If one forgets this effect of the mean field, Fermi acceleration 
at shocks looses rapidly its efficiency because the acceleration time 
increases like e 2 . One might think that the mean field increase the 
acceleration efficiency by maintaining particles in the shock vicin- 
ity during a longer time, as in perpendicular non relativistic shock 
acceleration. This is not so however, because the large advection 
speed of the downstream relative to the shock easily beats the trans- 
verse diffusion of particles, which can set in only after a time v" 1 . 



3 TRANSPORT OF SUPRA-THERMAL PARTICLES IN 
THE WAVE FRAME 

In this Section, all the quantities are measured in the frame of mag- 
netic disturbances, unless otherwise stated. In this frame the total 
energy of any particle is a constant of motion due to time translation 
invariance. 

In the upstream frame, the wave-numbers of Weibel modes 
obey k ± /k^ » 1 and %a) <k k ± , so that in the wave frame VJk!^ » 1 
as well. Regarding OTSI modes, one finds k ± ~ lq in the upstream 
frame, but Lorentz boosting to the wave frame gives predominance 
to the transverse wave-numbers. If, in a first approximation, the 
spatial dependence along the shock normal is disregarded, the nor- 
mal component of the generalized momentum also becomes a con- 
stant of motion: py + eA-(x,y) = C s1 . Thus the momentum of the 
particle is confined in a subset of the energy surface determined by 
Apu = eAAj, where AA. is the rms variation of the normal vec- 
tor potential. In short, the assumption of translational invariance 
along the Weibel filaments leads to an inhibition of momentum 
diffusion, see also I Jones et al. I dl998h for similar issues. This is a 
crucial point which directly impacts the efficiency of acceleration, 
which requires transverse scattering in the absence of a mean field. 
As we discuss further below, the transverse momentum is subject 
to large angular variations in the transverse plane, which leads to 
spatial diffusion. In order to obtain pitch angle diffusion, it is thus 
necessary to consider the full 3D dependence of the magnetic fluc- 
tuations; this scattering will lead to spatial diffusion along the three 
directions. 

Weibel modes propagate with a sub-luminal velocity v m along 
the beam (%(!> = k^v m ) and thus have a rest-frame in which parti- 
cles experience a static electric field and a static magnetic field. 
Similarly OTSI modes have a single sub-luminal phase velocity be- 
cause the wave number along the beam is selected by a narrow res- 
onance. For both cases the analysis of particle dynamics is suitable 
in the wave-frame, because this is the frame where the transport 



6 /. Plotnikov, G. Pelletier and M. Lemoine 



coefficients can be properly defined, where the distribution func- 
tion tends to become more or less isotropized, and fundamentally, 
this is the proper frame of scattering centers involved in the Fermi 
process. 

The Lorentz force acting in the wave-frame depends on the 
fields obtained by the Lorentz transform from the upstream plasma 
frame to the wave-frame: 

E\ = fin 

E' ± = y m (E ± +B m xB ± ) 
B\ = By 

B' ± = y m (B ± -B m xE ± ) 

In the case of Weibel modes, \E\ <k |6|, and for y m large enough, at 
first approximation, B' ± ^ y m B ± and E' ± = y„J3 m x B ± ^ B m x B' ± . 

In the case of OTSI modes, in the upstream frame the modes 
are almost electrostatic \B\ <k \E\, and E± ~ E\\ (oblique modes). 
Thus we neglected the contribution of Bp and B ± . The system is 
quite similar to the system derived for Weibel modes in their proper 
frame, since E' ± = y m E ± and B' ± m -y m B m x E±, which leads to 
B' ± a -B m xE' ± . We will analyze the particle dynamics first in a 2D 
approximation, meaning fijj — > 0, Bjj — > 0, and then in the complete 
3D configuration (fi,j = E' ± /y m ). 

When going to the wave frame, the perpendicular coherence 
length remains unchanged, tf ± = £ ± , while £ = y m C\\. This sec- 
tion considers quantities in the wave frame, as made explicit by the 
primes. 

3.1 Transport in 2D approximation 

The particle dynamics can be further simplified, when one neglects 
the longitudinal components (z direction) fijj et B'^ in a first approx- 
imation. Then the system reduces to: 



dp, 

= qE' x (\ + B m B z ) 



dPy 

dt 

dp, 



qE' y {\+B a fi z ) 



(21) 
(22) 
(23) 



= -qB m (B x E' x +ByE' y ) 

where the relation B' ± = —B m x E' ± has been inserted. Because 
K = ym(k z — B m u)lc) =i k-/2y m , for large y m the z dependence of 
the field can be neglected in a first approximation. The system is 
the same for both Weibel and OTSI modes. 

Because the turbulent fields depend on two spatial variables 
(x, y), there are two invariants: the total particle energy, given by 
the Hamilton function 



H = e(p) + q®(x,y) 

and the generalized momentum 

n z = Pw + qA z (x,y)/c . 



(24) 



(25) 



These two potentials are related to one another: A-(x, y) = <t>(x,y). 
Since the potential has a zero average and a finite rms value AC>, the 
particle proper energy e and its z-momentum py have well defined 
rms variations under the ergodic assumption: Ae = Ap\\c = eA</>. 
These relations determine confinement regions in phase space and 
forbid some diffusion processes. The norm of the transverse mo- 
mentum also has bounded variations. Indeed, considering the vari- 
ation from initial values, Se = e-eo,6p\\ = P\\-p\\,o,6p^ = p ± -p ± o, 
one exactly finds 



where we took into account the two previous invariants that imply 
Se = cSp\\. Now because p ± .o is small, this constraint allows large 
variations of the polar angle of p ± . Thus the variations of the energy 
and momenta are bounded, except for the angle of the transverse 
momentum that can vary randomly over the interval (0, 2n); those 
erratic variations of the angle can occur with the B z contribution to 
the transverse equations of motions that opens the phase space with 
another degree of freedom. This will have consequences on spatial 
diffusion, as will be seen further on. 

The transport properties are to depend on the particle momen- 
tum, and one may define a reference energy 



efi'f, 



(27) 



As discussed before, ajy' <sc 1 for supra-thermal particles, meaning 
that <k 6. Since these particles are strongly beamed toward the 
longitudinal direction, e - P||,o c initially, and p±$/P\\,o S; l/r s . 
In this limit, the transport of the supra-thermal particle can be 
described in configuration space as the random walk of a non- 
relativistic particle in the transverse plane, coupled to ballistic mo- 
tion along the longitudinal direction. In the transverse plane, the 
particle is described as non-relativistic because its transverse ve- 
locity v ± - p±/p\\ and, given that <K P\\$c, the particle cannot 
exchange a large fraction of its parallel momentum with transverse 
momentum due to the invariance properties in this 2D approxima- 
tion. 

The transverse motion, devoid of linear resonance but gov- 
erned by a continuum of Fourier modes, is thus characterized by a 
single nonlinear time 



<m = — 



tx 1 2P][,oc 



i /- 



(28) 



This is the time needed to cross a coherence cell for a particle that 
gets accelerated transversely in the transverse electric field. Over 
a coherence length, the electric field can be considered as constant 
and the particle receives a transverse kick cAp ± ~ e t . 

If the initial transverse momentum p ±0 <k e*/c, the spatial 
transverse motion can be approximated by 

1 eE ± 



Ax 



2 Pn.o 



ct 2 + u ± t . 



(29) 



since the kick remains much smaller than (e*p||,o) , v ±,o — 
p ±0 /p\\,o- The nonlinear time so defined is the time beyond which 
the nonlinear dynamics de-correlates the trajectories. It can also be 
considered as the time step for a random deflection of angle 0, in 
the transverse plane since each crossing of a coherence cell in the 
transverse direction is associated to a large variation of when 
p ±0 <k e*/c. Thus after n steps of size £ x , the trajectory has dif- 
fused such that 



<Ax 2 ± > - X -l\n , 



(30) 



with n ^ At/trn, which leads to a transverse spatial diffusion coeffi- 
cient 



4 ta 



Hl,o 



(3D 



2(eo - P\\,oc)6e = c 2 (6pj 2 



(26) 



where the numerical prefactor is subject to uncertainty. 

If P_l,oC » e*, the transverse velocity undergoes small varia- 
tion of its modulus in the crossing of a coherence cell, but it can 
undergo significant angle variations. This transverse quasi scatter- 
ing can be analyzed in two regimes: (a) when p ± ,oc <k (e*.P||,o) 
and (b) when p±^c s> (e*P\},o) but still p ± <k py . In the former 
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limit, the particle crosses a transverse coherence cell in a nonlin- 
ear time scale as previously, whereas in the latter case the crossing 
occurs on a linear timescale Let us estimate the transverse 

diffusion coefficient for both cases. In case (a), when crossing a 
coherence cell, the particle undergoes a small transverse deflection 
of angle (9, ~ e+lp L c. The scattering time is therefore longer than 
fa with r scatt ~ 6j 2 fa and diffusion in the transverse plane during 
Af > fscatt is such that 



At 

^scatt 



(32) 



This indicates that, in this regime of small deflection, a particle 
travels over a much larger distance than £ ± during a scattering time, 
namely (fscattAniKx- The transverse diffusion coefficient is thus 

1 l\ lp L c\ 2 

D ^ii\-)- (33) 

In case (b), the small deflection that occurs during the linear cor- 
relation time leads to a scattering time r sciitt ~ 9~ j ~ 2 C ± /v ± . And the 
transverse diffusion coefficient becomes: 



Pj_£ 



(34) 



3.2 Numerical results in ID-approximation 

The system of equations of motions in the fields is solved by the 
Bulirsch-Stoer algorithm l lPress etalJl 19861) , together with statis- 
tics over a sample of random phases and polarization directions 
of plane waves. The field E' ± is decomposed in plane waves with 
the constraint k' || E' ± , which implies that k' is in the plane (x,y), 
choosing as before the normal direction along the z-direction. In 
this subsection we set k'. = and E' z = in the wave frame. Sev- 
eral hundreds of modes are used in the setting construction. The 
system converges quite rapidly with ~ 10 3 particles. For simplicity 
we investigate only the case were the initial particle momentum is 
oriented along the parallel direction, p = p\\,oe z . 

All simulations use the following units: the spatial length unit 
is mc 2 /(eE' ims ) (E'^ is the root mean square of the electric field 
strength), the time unit is mc/(eE' rmi ), and momenta are expressed 
in units of mc. The coherence length £ ± ~ 1 in these units. 

Numerical results are presented in Figures 1 1 121 Fig. [T] de- 
picts the time evolution of spatial diffusion coefficients for different 
P\\f)l{mc) in a 2D approximate geometry. Green curves depict the 
time evolution in the parallel direction and blue curves represent the 
same quantity in the transverse direction. Displacements in the par- 
allel direction are clearly ballistic: (Ajc?) = uj t 2 - In the transverse 
plane, at least two regimes are identified: for t < fa, (Ax 2 ) oc f 4 as 
expected from particle acceleration by nearly constant electric field 
over its coherence cell, see Eq. 1291 for t > fa, (Ax 2 ± )/At reaches 
a plateau (i.e. diffusion) and its value has a power-law dependence 
on P|| i0 with slope of -1/2 (see the sub-panel of the figure), which 
fits well the result of the previous Section, see Eq.[3T] This behav- 
ior is a direct consequence of particle trapping in every coherence 
cell it encounters, with a "waiting time" per cell being equal to fa. 
The case where p ± is different from zero was also investigated. In 
this case there is no more trapping in the field coherence cell but 
random motions in the transverse plane. We obtain a diffusive be- 
havior on longer time-scales as discussed in the previous Section 
(e.g Eqs.[H[33}. 

Figure|2]presents the evolution of (pjf)/pjf„ (green curves) and 
(P±)/Pyo (blue curves) as function of time for different pp/Owc) 6 




Figure 1. Time evolution of spatial diffusion coefficients for different 
PWfiKmc) in 2D fields geometry approximation. Green curves: parallel di- 
rection; Blue curves : transverse direction [(x, y) plane]. Displacements in 
the parallel direction are clearly ballistic: (Ax 2 ) = v 2 Q t 2 . In the transverse 
plane, at least two regimes are identified: when t < fa, (Ax^ ) oc ( 4 as ex- 
pected from particle acceleration by nearly constant electric field over its 
coherence cell; and when t > fa, (Ax 2 ± )/At reaches a plateau (i.e. diffusion) 
and its value has a power-law dependence on p\\ q with a slope of - 1 /2 (see 
subpanel). Subpanel: Transverse diffusion coefficient (plateau) dependence 
on initial particle momentum p\\ rj. Dashed line follows the power-law slope 
-1/2. 



2D: (p. ) /p~ as function of time. Blue : ± dir. Green : II dir. 




Figure 2. Time evolution of {p 2 )lp 2 a for different values of p\\$l(mc), 
in 2D approximation. Green curves: parallel (z) direction; Blue curves : 
transverse direction. Subpanel: Asymptotic values of (p 2 )/p 2 as function 
of p\iQ. (Py)/Py ' s slightly inferior to 1 independently from py o (green 
squares). In transverse direction (blue circles) (p^/p 2 oc p^. 



[1, 10 4 ]. There is no diffusion in momentum space with parallel 
component being bounded by the n- invariant, while the transverse 
one varies as (A/>5.)/(mc) 2 oc p^ . As generalized momentum is 
constrained by this n L invariant, we expect a transverse momentum 
gain (6p 2 ± ) = eE'J ± p l0 /c. 

Finally, we can mention that the energy gain (Ay) is indepen- 
dent of pio and corresponds to the amount of energy brought by the 
rms electric field potential. It can be expressed as (Ay) = eE' ± t ± . 
The energy variation from e(/>n,o) to e(pi )+qE'f ± takes place when 
< t < fa. For t > fa, the particle energy remains constant. 
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3.3 Transport in 3D-fieIds 

We consider now the particle transport over a longer time with mo- 
tions in the z direction, which involves the single resonant mode 
with k z = 6~ l in the case of OTSI and a continuum of small wave 
numbers in the case of Weibel modes. In the wave frame, the lon- 
gitudinal coherence length = y m C\\, and l\\ ~ l x for OTSI modes 
in the background plasma frame, so that (',, y m i' x - For Weibel 
modes, £\\ > i x and therefore l'Jl' x > y m - 

The invariance of the generalized z-momentum is now bro- 
ken and the particle momentum can diffuse in all directions. One 
limitation is the phase-space confinement due to the total energy 
conservation, that forbids energy diffusion. Another is that the lon- 
gitudinal coherence length is significantly larger than the coher- 
ence length in the transverse direction. With respect to the above 
2D analysis, we thus expect a change of regime when the motion 
in the z direction is felt by the particle, which occurs on time scales 
much larger than fjj/c, when a particle has crossed many coherence 
cells along the z-direction. The randomization of the longitudinal 
component of the momentum is expected over a time t z , with 

t „ f '' f gggf f „ { \\ £ '± (Plivrf (35) 
c [eE'J'J t' x c \ e* ) 

For both OTSI and Weibel modes, E' z i!, - E' x l' x because in the wave 
frame, V x E' = 0. The above is a linear estimate of the scattering 
timescale in the longitudinal direction, however simulations at low 
energies suggest that the dependence on ft. differs from this latter 
We will take into account that lack of knowledge by introducing a 
factory measuring the delay, compared to the 2D correlation time, 
in the full development of the 3D-dynamics. It will turn out to be 
important for the performance of the acceleration processes. 

Over intermediate timescales, we expect to recover the previ- 
ous results about the transverse diffusion and no diffusion in the 
longitudinal direction. But on time scale longer than one expects 
a 3D diffusion with 

D\\ = fD x . (36) 
x 

As before we assumed that the deflection over a coherence cell is 
small with S6 ~ e*/pu ,oC. The description of the expected behaviors 
of these diffusion regimes is provided in the figures shown in the 
next paragraph. 

3.4 Numerical results in 3D-flelds 

We take the same configuration as in the 2D case but with a finite 
value of E'., = E' x /y m in the case of OTSI, with k'^/k' ± = l/y m . To 
simulate a Weibel turbulence, we include an ensemble of modes, 
consistent with the original wave equations of Weibel modes, 

E;^) = --£(k' x xB' x )-f) m , (37) 
k x 

which insures the conservation of the total energy of each particle. 
For this ensemble of Weibel modes, kyk' ± takes small values up to 
some parameter k < 1; consequently, E'JE' X ~ k. Most numerical 
calculations were done with the same time scales as in 2D case. 

Figure [3] depicts the time evolution of spatial diffusion coeffi- 
cients for different p^ /(mc) in 3D OTSI fields, assuming y ra = 30. 
Green curves correspond to the parallel direction and blue curves 



to the transverse direction. Displacements in parallel direction re- 
main ballistic on numerical time scales and in transverse direction 
2D-like diffusive behavior (see Fig. |T} disappears gradually when 
f„i < t < t z . It will be seen further that the diffusive behavior is re- 
covered in all directions on longer time scales (i.e. t » t : ). Figure|4] 
presents the evolution of (p?}/p? in time for p^ /(mc) e [1, 10 4 ]. 
Comparing to Fig. [2] one observes that {p\)l p\ a keeps the same be- 
havior as in 2D and the transverse components begin to rise at later 
times; longer time simulations are needed to explore its asymptotic 
behavior. Interestingly, a different behavior is observed for parti- 
cles with momenta satisfying f„i > t\\ fi : in this 3D OTSI turbulence, 
electromagnetic fields reverse every half coherence length along 
the longitudinal direction, due to the periodicity in that direction; 
therefore particles with t n \ > fy L . execute oscillations in the trans- 
verse direction, with ballisitic motion along the longitudinal direc- 
tion. Such particles are then confined in the transverse plane be- 
cause their motion is reversed before they have time to experience 
a decorrelated field in the transverse plane. This is particularly im- 
portant with respect to acceleration efficiency, since such particles 
would not return on a short timescale (in the absence of a mean 
field, of course). One can check that the condition /„i > fp iC amounts 
to pn,o > e*y m /2 for p ± o = 0, which is easily satisfied, or in terms 
of initial transverse velocity, Pj.,o/Ph,o < 2/y m , which is also generi- 
cally satisfied for the supra-thermal particles. This means that OTSI 
turbulence is inefficient from the point of view of scattering supra- 
thermal particles away from the longitudinal direction. 

In order to make an assessment of our estimates from sub- 
section [33] a simulation with enhanced integration time was per- 
formed in the case where p^/mc = 1. The result is presented in 
Fig. [5] where different statistical quantities are plotted as a func- 
tion of time. A slightly smaller value y m = 10 is adopted here 
to reduce the characteristic t- time. On time scales larger than t- 
the spatial diffusion is recovered in all directions, all momenta are 
isotropized and energy gain is equal to the field rms energy. Note 
that the ratio between spatial diffusion coefficients in the parallel 
and transverse directions is not equal to y m as expected in Eq. [35] 
its value is 55 > y m and scales as 0.6y m when y m is varied explic- 
itly in our simulations. It remains uncertain if this scaling depends 
on the choice of too low energy particles (pio/(mc) = 1) and if it 
remains the same for highest energies. Direct simulations for the 
former case, with large enought integration time, are numerically 
prohibive and are subject to severe numerical errors. In the follow- 
ing, we encode this uncertainty in a parameter^, such that^- = y m 
if the linear value given by Eq. [35] were to apply, but x ~ 7m as 
indicated by the simulations. 

Finally, in the Fig.[6]we present also the same simulation but in 
the case of 3D Weibel-type fields with k' x /k',, = 10, corresponding 
to k = 0.1. The general behaviour is similar to that observed for 
transport in OTSI turbulence. In particular, one recovers a scaling 
t- oc k~ 2 and D\\/D x oc k~ 2 instead of the linear estimates given 
in Eqs. 1351 and 1361 which suggest a scaling in k~' . As before, we 
encode this uncertainty with a factory, so that D^/D ± m x< m e 

simulations indicating^ = (^y/^l) • 

As expected, the energy does not undergo a diffusive behavior. 
After a limited gain, the particles keep a constant momentum. This 
is in agreement with theoretical predictions: stochastic acceleration 
is not seen, transverse heating is bounded, the energetic particle 
beam is broadened in the transverse direction, but the distribution 
remains anisotropic. 




Figure 3. Transverse spatial diffusion in 3D modes of OTSI type. This 
graph depicts the time evolution of spatial diffusion coefficients for differ- 
ent pi\ o/(mc) in 3D fields. Green curves: parallel (z) direction; Blue curves: 
transverse direction [(x,y) plane]. Dashed vertical line delimits the linear 
coherence time in z direction f|| c = y m /^/c. Displacements in the paral- 
lel direction remain ballistic on numerical time scales, since the integra- 
tion time is too short to probe the diffusive behavior in this direction. In 
the transverse direction, 2D-like diffusive behavior (see Fig. (T)) disappears 
gradually when t B \ <t <t z . 



Figure 5. For OTSI modes, time evolution of different quantities in the 
case pno = mc with enhanced integration time. Solid curves show the av- 
erage (p?)/p^q, red color for parallel direction and orange color for the 
transverse one. Dashed curves show spatial transport coefficients (Axj)/At. 
Green color for parallel direction and blue for transverse direction. Dot- 
dashed black curve shows (Ay). Vertical dashed lines indicate tree carac- 
teristic times relevant for particle dynamics. As expected, spatial diffusion 
is present in all directions on time scales much longer than t-. Above the 
diffusion time all momenta reach isotropy: (p 2 x ) = (py) = {p}) ^ p^ „/3. 



3D : ( p )'p ||0 iL5 function of time. Blue : ± direction. Green : II direction 
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t 

Figure 4. Time evolution of (p^)/p^ for different values of py o/(mc), in 
3D-fields. Green curves: parallel (z) direction; Blue curves: transverse di- 
rection. Dashed vertical line delimits the linear coherence time in z direction 
t\l,c = yml±/c. Comparing to the Fig. (Pp/pf\o keeps the same behavior 
as in 2D and the transverse components begin to rise at later times but longer 
time simulations are needed to explore its asymptotic behavior. A different 
behavior are observed for particles with momenta satisfying /„i > ry c , which 
are confined in a coherence cell in the transverse plane, see the text for de- 
tails. 
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Figure 6. Same as in Fig. \5\ but for Weibel 3D modes. Here, we take 
kt/k? ± = 1/10. This choice is dictated by numerical time limitation in or- 
der to observe the development of a 3D regime. As expected, spatial dif- 
fusion is also present in all directions on time scales much longer than 
t z . Above the diffusion time all momentum components reach isotropy: 
(p\) = (pI) = (p\) — p\ o/3- We remark that the evolution is very similar to 
that obtained with OTSI modes, despite the finite range oik', wavenumbers. 



4 ELECTRON HEATING 

As we discussed in Section[2l the wiggler parameter for the back- 
ground electrons in the proper frame of the micro-turbulent mode 
is very large, actually a » Y, with y' the typical Lorentz factor 
of the electrons in that frame. Since the modes carry a transverse 
electric field of the same order of the magnetic field, this offers a 
promising source of preheating in the shock precursor. 

We first consider the effect of the electric field of the Weibel 



waves in their proper frame. The electron temperature temporarily 
increases, as long as their energy is smaller than MC>' = eE'C ± , be- 
cause, in this frame, the total energy of each particle is conserved, 
as discussed before for supra-thermal particles. This limiting en- 
ergy is 

6 = eE'€ x = ^B 2 7m-^-m p c 2 . (38) 

Reverting to the background plasma frame, that energy corresponds 
to a temperature which is a sizeable fraction of the proton energy 
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,1/2 *i- 2 
ft Jm-^m v C , 



(39) 



assuming that most of the electron heating is distributed along the 
transverse direction. This transitory heating process is not in a dif- 
fusion regime, it is rather a direct linear acceleration in the coherent 
electric field of a coherent cell. Using Eq.[23] one can easily check 
that the energy e' is obtained over a typical linear timescale £±/c. 
This fast heating in Weibel waves is a particular case of a situation 
where thermal electrons undergo strong relativistic motions in the 
waves, reaching relativistic temperatures T e ~ am e c 2 . A similar 
process is at work in the OTSI turbulence. 

Because the transverse coherence scale is everywhere smaller 
than the precursor length scale, on which £b varies, the above 
fast heating process brings forward a picture in which electrons 
are nearly instantaneously heated to the local temperature given 
by Eq. [39] above, scaling as ^ 2 . Of course, as the electrons near 
equipartition with the ions, one may expect the above heating pro- 
cess to saturate, notably because the oscillation parameter then be- 
comes smaller compared to y', as discussed in Sec. [2] Far from the 
shock, f B « 1 hence the electrons are heated to temperatures well 
below m p /m e (upstream frame), while closer to the shock, the tem- 
perature rises. Interestingly, the simulations of Sironi & Spitkovsky 
(2011a) indeed show a gradual evolution of the electron tempera- 
ture over the length scale of the precursor. The above preheating 
process provides a concrete physical mechanism for this picture. 

The previous estimate of the electron temperature allows to 
get a better knowledge of the filament size. The filament transverse 
equilibrium can be obtained from a simple pressure equilibrium 
argument, according to which: 5nT e + SB 2 /8/r = 0. This leads to 

\6n\ T e 



ft 



(40) 



and the existence of filaments implies a significant density trough 
such that \6n\ln ~ 1. Since 



e = ft 7m— m p c 

0\ 



(41) 



the coherence length and thus the filament size is such that £ x ~ 
y ra '£g /2 (5j and the temperature 



ftm p c 



(42) 



thus y e fi ~ £b and the electron inertial scale (accounting for tem- 
perature) compares well with the coherence or filament radius: 

6 e ~ <f B /2 <5. • (43) 
The filament motion in the hottest phase is then determined through 

7m ~ (fti-ft) 1 ' 4 for close equipartition 

~ Q l/2 for ft < f n , (44) 

which should thus increase with the distance to the shock. 

Close to the shock front, the filament scale A in the normal 
direction to get reflection is 

1 



£ m y 

?B Ym 



Si 



(45) 



so that 




A ~ 


(1 






or 




A ~ 


(1 







L/4 



1/2 



Si , for ft < ft 



(46) 



If we furthermore postulate that A corresponds to the efolding 
length scale ( g = ft- 1/2 (5i, then we find a simple relation between 
the two fundamental conversion factors: ft ~ ft in the first case, 
ft ~ ft. in the second ca se. However the numerical works by 
ISironi & Spitkovskvl j2oTT3) indicate a better performance of mag- 
netic field generation with £ B ~ ft in the case of a proton electron 
plasma; which suggests a steeper barrier: A ~ S\. Only the cosmic 
rays fraction ft r remains undetermined. 

For background electrons this heating occurs in a very short 
time t ± /c, a few ftV 2 w pi ', which is shorter than the growth time. 
The electron pre-heating has an important feedback on the insta- 
bility. Although it does not saturate the instability, it reduces the 
phase velocity of the Weibel modes, at least because it increases 
the small mass ratio and determine the condition for the reflection 
of incoming protons. 

Considering now OTSI modes, we have argued in section 4 
that the heating process is similar, and we expect that it stops with 
the saturation of the OTSI when the electrons achieve a relativistic 
temperature. 

So far, we have discussed the case of a nearly unmagne- 
tized shock. One may recall that at high magne t ization, the wake 
field develops in a p + e~-pla sma dHoshinoll2008l . lLvubarsky||2006l 
ISironi & Spitkovsky||201 13) . In the frame of the wake wave, there 
exists also an energy invariant that leads to a relativistic tempera- 
ture for the electrons, but not to acceleration. 

Across the shock front, the proton heating follows from the 
mixing of the different proton flows carrying energy of order 
T s m p c 2 in front frame. More puzzling is the issue of electron heat- 
ing. In the previous paragraph we saw that a significant preheat- 
ing of the electron is expected. The final stage of electron heating 
across the shock front is necessarily related to an effective loule 
effect due to scattering. The scattering frequency of thermal parti- 
cles in the shocked flow is v s , ~ ft<w p i and the magnetic diffusivity 
v m = rjc 2 /(4n) = v„5 2 ~ ftcft. The typical length of Joules dissi- 
pation is thus tj = 3v m /c ~ 3ft(5j, which is a quite short distance 
for particle thermalization. A more detailed estimate is obtained 
by looking at the absorption of each Fourier mode, which leads to 
£j(A) — 3/l 2 /ft<5j. All the magnetic energy that has been generated 
in the precursor is thus dissipated in the electron heating, the elec- 
tron temperature remains in sub-equipartition with that of protons, 
T e ~ £ B T s m p c 2 . It turns out that both processes, pre-heating and 
Joules heating, have a similar contribution to the electron tempera- 
ture. 

Finally we have the estimates of all the scales involved in this 
physics, measured in the ambient frame they are: 

• The Weibel instability growth length: <5j/ftV 2 

• The background electron heating scale: ft^ 2 (5i 

• The transverse scale of maximum instability in the hot phase: 

S e ~ 

• The transverse coherence length C ± , also the filament radius: 

e± ~ £ 2 Si 

• The filament length: A ~ 6-J^'J 2 if ft ~ ft\; A ~ 6, with 
ft ~ ftr- 



5 APPLICATION TO ACCELERATION AT A 
GAMMA-RAY BURST EXTERNAL SHOCK 

Section[3]provides the tools required to discuss the residence time 
of particles upstream of a relativistic shock. As discussed in Pel- 
letier et al. (2009), one must compare the diffusion timescale in 
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the micro-turbulence with the timescale associated to rotation in 
the background field and keep the shorter of the two. Let us discuss 
here the implications of the micro-turbulence. The Weibel filaments 
are likely the best sites of scattering, since the OTSI modes, at least 
in their linear description, lead to confinement of high energy par- 
ticles along the shock normal. 

The residence time of supra-thermal particles returning from 
the shocked plasma in the upstream flow is that corresponding to a 
deflection of an angle 1 /T s after which the shock front catches it up 
again, and this fixes the acceleration time of the Fermi process. This 
has to be estimated in the filament frame first, in which the angle 
is l/7f| m , with 7f| m the relative Lorentz factor between the shock 
front and the wave frame, as determined by Eq. [9] The influence 
of a possible background magnetic field on the return timescale is 
discussed further below. For a particle of energy e ra in this wave 
frame, the residence time is 



■x- 



1 

f|m 



(47) 



where ^ is a factor large compared to unity. It accounts for the 
fact that the decorrelation time in the longitudinal direction is much 
larger than that in the transverse direction, as discussed above. For 
OTSI modes, the simulations indicate x ~ 7 ra > although the con- 
finement in the transverse plane leads to very ineffective scatter- 
ing; scattering is rather provided by the Weibel modes, for which 
X ~ y 2 a k±/k\] > y 2 v in terms of the wavenumbers of the instability 
measured in the background plasma rest frame. 

Going to the front frame, one finds a residence timescale 



fres.f 



i 

1 fjm 



(48) 



with = (y m /T s )eB v £ ± in terms of the turbulent magnetic field in 



rf[m(l ~A|m) e m- 



•A 



(49) 



c \eB\fi ±/ 

The fact that the scattering takes place in a frame moving at 
high speed shortens the residence time, but this gain is mitigated by 
the anisotropy of the turbulent modes, which induces the^ factor. 

5.1 Electron acceleration at relativistic shocks and radiation 
(GRBs) 

We have developed arguments in favor of an efficient heating of 
the electron fluid by the micro-turbulence, which confirms the idea 
that the electrons could likely reach a sub-equipartition temperature 
at relativistic shocks. For instance in GRB termination shocks that 
give rise to the afterglow emission, they could achieve a tempera- 
ture of a few tens of GeV. Indeed the proton temperature is very 
high at the beginning of the afterglow, and we have 

r s 

3V2 

which corresponds to a few tens of GeV. Intense short scale mag- 
netic turbulence develops because the interstellar magnetization pa- 
rameter is very low, cr ~ 10~ 9 . 

What kind of radiation can be expected in such small scale 
field, much more intense than the mean field? This depends on the 
wiggler parameter a, now measured in the downstream frame: 



T, < T„ 



(50) 



eBic ,i/2 r "2L 
m e c l m e 



(51) 



This parameter measures the capability of the magnetic force to de- 
viate a relativistic electron of Lorentz factor 7 by an angle 1 jy (this 
is the reason for which y does not appear in the definition). When 
a > 1 the magnetic field produces a single deviation of the electron 
in the emission cone of half angle I/7, whereas when a < 1 the 
electron can undergo several wiggles in the emission cone. When 
a is large, the emission behaves like a normal synchrotron radi- 
ation in a mean field, except that there is no polarization . When 
a is small, the emission is of jitter type dMedvedevll2000l) . In the 
present case, the large wiggler parameter ensures that the emis- 
sion caused by shocked and accelerated electrons at a relativistic 
shock is synchrotron-like; the analysis of the emitted spectrum may 
provide a diagnosis of the magnetic turbul ence although the de- 



partures are expected to be moderate (e.g . Kirk & Revillel 201ol; 
Reville&KirklhoiOl ; iFleishman & Urtievl 120101 :1 Medvedev et al.1 



201 ll) and likely dominated by the decay dynamics of the mi- 
croturbulence downstream, which imply that particles of different 
Lorentz factors cool in regions of different magnetic field strengths 
dLemoinell2012l) . 

As for the supra-thermal electrons, we find a similar estimate 
of the maximum Lorentz factor, measured at shock front, achieved 
agains t synchrotron loss than the one derived by iKirk & Revillel 
d2010h when y m ~ T s , except of the dependence in y m : 



4?rg 2 f c 7m 
ajm e c 2 x^s 



7 x 10 6 



*r 8 



1/3 



1/3 



/ i\- 1/6 / >" m 



(52) 



To obtain this result we took into account the electron wandering 
upstream where it experiences a level of turbulence comparable to 
the downstream one, as measured in the front frame, because the 
extension of its trajectory upstream is much shorter than the high 
energy proton trajectories that shape the precursor. The correspond- 
ing maximum energy for the emitted photons received by an ob- 
server is 



1-4/3 - 



2 > 3 v 2 ' 3 m c 2 

7m '"pi- 



(jmr 3 .)- 1 / 6 a f 



3 x ^B,-2r s ,2.5"C 



1/2 7, 



2/3 



GeV. 



(53) 



with the usual notation £ B _ 2 = £b/0.01, r s2 .5 = r s /300 and 
no = rc/1 cirT 3 . This relatively high performance is due to the pro- 
ton dynamics; otherwise with electron dynamics, only the MeV 
range would be reached. 

The above estimate of 7„ lax balances the acceleration timescale 
against the timescale for synchrotron losses in the turbulent field. In 
principle, one should also include inverse Compton losses, which 
puts strict const raints on the return tim e scale in the upstre am frame, 
as discussed bv lLi & Waxmarj ( l2006t) . lLi & Zhao! fcoi ll) . This re- 
quires to use the scattering timescale discussed above and follow 
the proper treatment of Klein-Nishina suppression given in these 
studies; this task is left for future work. 

The above discussion considers an unmagnetized shock. In 
the presence of a background magnetic field, return in the up- 
stream can be achieved through the rotation by an angle 1/T S in 
the background field. The return timescale then corresponds to 
fres.oif ~ e[f/(r s eBo) as measured in the shock front frame, so that 
at the maximal Lorentz factor determined by Eg. 1521 



'res,0|f 
fres|f 



2 7» 



2/3 



£s,-2Bo,-5' ! 



2/3 r l/3 



s,2.5 



(54) 
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with 60,-5 = Bo/10/jG. This implies that the background magnetic 
field starts to dominate the dynamics of the highest energy electrons 
as soon as the ratio y m /x becomes significantly smaller than unity, 
for the above fiducial values. Among others, this guarantees that 
GeV photons can be produced, independently of y m /x- 

Thus a single synchrotron-like spectrum extending up to sev- 
eral GeV, even possibly a few tens, can be expected and thus is 
compatible with observations. From that point of view, the perfor- 
mance of relativistic shocks for electron acceleration and radiation 
can be regarded as quite high. The conversion factor into radiation 
is f ra d ~ ^BO"T«''s<7e> (r s the external shock radius), and at the be- 
ginning of the afterglow £ rad ~ £ B ~ 10%. 



5.2 Relativistic shock and supra-thermal protons in GRBs 

As mentioned previously, because the scattering time and thus the 
acceleration time increase with e 2 , the Fermi process at relativis- 
tic shocks is not expected to be an efficient accelerator of protons 
towards the highest energies. That acceleration is limited by expan- 
sion losses of time scale r s /c in observer frame. At the termination 
shock of GRBs at the beginning of the afterglow, the maximum 
energy achieved when the residence time upstream balances the 
expansion time, is 



2zr s 2^ 



1/2 



x 10 15 xZl 7, 

X 



I -:<r., 11,, eV 



(55) 



Again, the above holds for an unmagnetized shock. The perfor- 
mance can be improved if one takes into account a background 
magnetic field, which leads to regular rotation and a shorter return 
timescale at these energies. The maximum energy can then be writ- 
ten £ max = ZT s eB r s ~ 0.3 x 10 7 x ZGeV. Thus although an 
energy of order 10 16 eV is achieved, the result is far from the range 
of so-called ultra-high energy cosmic rays. 



Here we assumed that r s a 1 so that B 2 /4n a £bPoc 2 . Therefore 
it is quite valuable to consider mildly relativistic shocks in ultra- 
relativistic flows. 

The jets of FR2 radiogalaxies have a power P iet = 10 44 - 
10 46 er g/s, hence they might produce UHECRs jPtuskin et al.l 
fcOllh ). GRBs flows carrying internal shocks can be even more 
efficient accelerators, depending of the duration of the flow for a 
given energy. For instance, a GRB of apparent isotropic luminos- 
ity Pj c t/(f 2 /4) ~ 10 52 erg/s (9j being the half opening angle of the 
flow > l/r ; ) with T ; ~ 100 may produce particles with energy as 
high as Z x 10 21 eV for a similar conversion factor £b = 0.01. Thus 
mildly or sub-relativistic shocks in a relativistic flow are more ef- 
ficient accelerators of protons than ultra-relativistic shocks and are 
excellent candidates for being sources of UHECRs, thanks to the 
magnetic field amplification at shocks. 

Now if we consider FRl-jets, they are less powerful than the 
FR2, with P Jct = 10 42 - 10 43 erg/s and they are sub-relativistic. The 
maximum energy achieved in those jets is at least one order of mag- 
nitude lower than that achieved in FR2-jets. In principle, this might 
be compensate d by the generation o f UHECRs in the form of iron 
nuclei (Z=26) l lAloisio at alj|201 lh . For instance, we have fairly 
precise VLBI measurement of the radius of Centaurus A at several 
distances from the nuclei, together with estimate of the magnetic 
field responsible for synchrotron emission; this leads to an estimate 
of the product Brj ^ 2 x 10 ~ 3 G pc. The Hillas criteriu m leads to 
Emax - Z x 2 x 10 18 eV (see lLemoine & Waxmar]|2009l) . which is 
consistent with a significant amplification of the magnetic field at 
internal shocks. The mos t recent VLBI map of Centaurus A, done 
by the TANAMI project dMiiller et"al]|201 ll) , reveals structures in 
the flow, suggesting internal shocks. The recent data analysis of 
Pierre Auger Observatory suggests an enrichment of the UHECR 
spectrum with heavy elements at higher energies and also a possi- 
ble correlation with Centaurus A. However it is very difficult to un- 
derstand the observed pattern of anisotropy i f one assumes that the 
obser ved particles are heavier than hydrogen dLemoine & Waxmanl 
120091) . In this sense, a contribution of protons coming from more 
remote and more powerful sources appears unavoidable. 



6 ABOUT THE ORIGIN OF UHE-COSMIC RAYS 

Precise performances of mildly or sub-relativistic shocks are not 
yet known and require more numerical simulations. However 
some reasonable guess are permitted by extrapolating what we 
know about the two extremes: non-relativistic and ultra-relativistic 
shocks. The main guess is that we can expect a magnetic field 
amplification at shocks wit h a convers ion factor f B = 1 - 10%, 
occurring in MHD regime jBeiill2004l) without severe limitation 
due to the super-luminal configuration, especially for oblique inter- 
nal shocks (termination shocks in the hot spots of FR2 jets might 
be super-luminal). These assumptions can be applied to internal 
shocks of AGN jets (in particular in Blazars jets), and to internal 
shocks of GRBs; the limitation of the acceleration of protons and 
nuclei is essentially the escape of particles at the edge of the flow. 
For a relativistic jet of radius rj, with a large bulk Lorentz factor 
Tj and a kinetic power Pj et = poc 3 r 2 7ir 2 (po being the co-moving 
density), the maximum energy given by assuming acceleration at 
the Bohm limit, i.e. an acceleration timescale f acc = r^/c in the co- 
moving frame, the maximum energy in the observer frame reads (in 
agreement with Lemoine & Waxman 2009): 

£ max - rjZeBrj - Z x 10 19 ( f*, f je ' ) eV . (56) 
\ 10 2 10 45 erg/s / 



7 CONCLUSIONS AND PROSPECTS 

The triangular dependence of collisionless shock structure with a 
reflecting barrier for a part of incoming particles, with generation of 
supra-thermal particles and the generation of magnetic turbulence 
is a successful paradigm that applies to astrophysical shocks, both 
non-relativistic and relativistic. Numerical and theoretical works 
are making significant progress for both understanding the physics 
and providing quantitative results useful for astrophysical investi- 
gations. This includes not only the spectrum index and cut off of the 
distribution of accelerated particles, but also the conversion factors 
into cosmic rays, magnetic turbulence and radiation. We are likely 
seeing only the beginning of such studies, that require more PIC 
simulations and new types of hybrid codes involving relativistic 
MHD coupled with PIC codes for cosmic rays. 

In this paper we have presented new theoretical investigations, 
namely the crucial effect of the motions of magnetic disturbances, 
the detailed analysis of the contribution of OTSI and the issue of 
the pre-heating of background electrons. The electron preheating is 
insured by a motional electric field in the wave-frame, the inten- 
sity of this field is almost the same of the magnetic field. Particles 
experience scattering and heating at the same rate. 

Despite that the Weibel instability generates magnetic fila- 
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ments whereas OTSI generates almost electrostatic waves in the 
background plasma, they behave similarly in their proper frame, in 
which they are composed of an electro-static field and a magneto- 
static field of almost the same amplitude. The supra-thermal par- 
ticles scatter in these fields while the background electrons are 
pre-heated up to sub-equipartition in such fields, because they ex- 
perience relativistic oscillations. The preheating is particularly ef- 
ficient: within a coherence length of the perturbations, it heats 
the electrons to ~ £ B m p c 2 , in which f B should be understood as 
the local (position dependent) fraction of energy density stored in 
the electromagnetic component. Because the coherence length is 
much shorter than the size of the precursor, this brings forward 
the picture in which the electrons are instantaneously heated to 
the above temperature, so that their temperature rises gradually to- 
wards near equipartition as they approach the shock front, a pic- 
ture which appears i n sati s factory agreement with the r esults of 
ISironi & Spitkovskvl J2009h . lsironi & Spitkovskvl feOl lab . 

Electron pre-heating modifies the generation of micro- 
turbulence: it saturates OTSI and slows down the propagation of 
Weibel modes. So we envisage that the nose of the precursor con- 
tains fast propagating Weibel modes and then, closer to the shock 
front, relativistic thermal electrons that quench OTSI, slow down 
Weibel modes, enlarge the characteristic scale. The reflection of 
incoming protons is made possible under that condition of hot 
electrons with an enhancement of magnetic energy density up to 
a fraction of pc 2 over a distance of a few proton inertial scales 
Si. Thus the OTSI instability, although it grows a little faster than 
the Weibel instability, plays a transitory role in the precursor of 
an ultra-relativistic shock; although it plays a role similar to the 
Weibel instability, it saturates when the electrons become relativis- 
tically hot. The OTSI modes are probably not transmitted down- 
stream, where they would propagate as electromagnetic waves of 
a quasi-vacuum. Nevertheless they are definitely active in the cold 
phase of the shock precursor, like Buneman instabilities that also 
preheat efficiently the electrons. 

The motions of magnetic filaments relate to the 3D description 
of the modes, and the 3D description is crucial for the development 
of the scattering process. Clearly, because of the filament nature 
of these magnetic structures, the time scale for a full development 
of 3D scattering is essential. Indeed the performances of the accel- 
eration process can be reduced significantly (as described by our 
factory) as t he time required to reach 3D-scattering increas es. The 
simulations dSpitkovskvll2008d.lsi"roni & Spitkovsky||2011al) of the 
Fermi process in 2D probably involves mirror effects on the shock 
front rather than actual upstream/downstream scattering, especially 
at the "low" energies corresponding to the first Fermi cycles probed 
by these simulations. To probe the 3D scattering regime discussed 
here, one would need 3D PIC simulations with very long integra- 
tion timescales, in order to accelerate particles to energies such that 
their Larmor radius in the turbulent field becomes larger than the 
coherence length. 

The combination of Weibel and OTSI turbulence can produce 
a kind of relativistic second order Fermi process upstream, in ad- 
dition to the effect of their own electric field. However the interac- 
tion time is limited by the residence time of supra-thermal particles 
upstream; therefore a significant improvement of the acceleration 
process is doubtful. 

Shocks in AGN, blazar jets, or in the internal flow of GRBs 
are mildly relativistic and therefore not subject to the severe re- 
striction imposed to the Fermi process by the mean field like 
what happens in the ultra-relativistic termination shocks of pul- 
sar wind nebulae and of gamma-ray bursts. In pulsar wind nebu- 



lae, reconnections likely contribute to injecting high energy par- 
ticles in the shock and a supra-thermal tail with a hard compo- 
nent may be generate d dLvuba rskv 20031: IPetri & LvubarskvlkoOTl : 
ISironi &"Sp itkovskv 2 011b ). In the termination shock of a GRBs, 
the magnetization is very weak and Fermi acceleration becomes 
operative when it decreases below some critical value, because 
the conditions for generati ng magnetic turbulence are fullfilled, a s 
suggested by observations dLi & Waxm an 2006; Li & Zhao: 20 lH ). 
When the termination shock propagates in a star wind whose den- 
sity decreases like 1/r 2 , the Fermi process may be, in a first step, 
quenched by the mean field, if the magnetization is large; then as 
the magnetizat ion decreases, it starts and generates very high en- 
ergy electrons dLemoine & Pelletierl201 lbl) . When the shock prop- 
agates in a normal interstellar medium of very low magnetization, 
the Fermi process works at the very beginning. Fermi process at an 
ultra-relativistic shock is a very efficient accelerator of electrons, 
the scattering rate at low energy being very high. And because it 
generates also a very intense magnetic field at short scales thanks 
to the incoming proton energy, it produces a very intense and ex- 
tended synchrotron-like spectrum up to GeV energies. 
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APPENDIX A: WEIBEL INSTABILITY WITH A 
NON- VANISHING PARALLEL WAVENUMBER 

Despite detailed analysis of the relativistic beam instability i n 
Weibel regime, like for in stanc e in Achterberg & Wie rsmal d2007l) . 
IShaisultanov et al.l l l201ll) . and lBret et al.l J2005t) . the motion of fil- 
aments has not been given attention so far. We will emphasize this 
issue for cold and relativistically hot electron fluid. 

For a cold background plasma pervaded by a cold beam, the 
wave system is described by the matrix: 



k 2 C 2 , 



(Al) 



where x\ are the components of the susceptibility tensor of the 
beam plasma given by 
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We chose a beam velocity direction along z-axis, v b = v b e. and a 
wave vector k = k x e x + k z e z . Electromagnetic waves polarized in 
the y-direction are decoupled. The wave system of interest reduces 
to second order, such that 
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And the dispersion relation is 

D(k, oj) = A XX A---Al z = 

In the case of Weibel instability, oj 2 <k a) 2 e , k 2 6 2 
k 2 8 2 <K 1, we find: 
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with k ± = k x . This is the result for a cold background plasma and a 
cold beam. When one takes into account a dispersi on of the beam 
withi n an angle l/r s , as previous ly investigated dRabinak et al.l 
1201 ll : iLemoine & Pelletien |2oTT3) . the instability is quenched 
when T s < (tjcyfiy 112 . However the instability is restored when the 
background electrons are sufficiently hot with both the same growth 
rate and the same frequency. In the case of OTSI, electron heating 
up to relativistic temperature tends to quench the instability because 
the modes become super-luminal and thus the resonant interaction 
is no more possible; this effect is however delayed by the fact that 
the frequency is negatively shifted by the beam, which lowers the 
phase velocity. 

Let us now study Weibel's instability properties for a rela- 
tivistic electron temperature. As the electrons can come close to 
equipartition, it is essential to account for the response of the back- 
ground ions. Let us write the components of the wave tensor: 
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k 2 V - > ) k 2 

where the dielectric coefficients for relativistically hot electrons and 
cold protons in the low frequency approximation are: 



%i 1 / n oj\ 

v -^ + ¥KX + l 2kc) 



or. 



,ii kc 

4 £J 



(A6) 
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When considering relativistically hot electrons, it is convenient to 
write the Debye length X De such that 

T 1 
^Dc - — e ~T = Ptf with p = -y e ii and y e = 1 + 3T e /m e c 2 . 
Ame 1 3 

Whereas the Landau contribution (imaginary part) is a small cor- 
rection in the longitudinal response, it is dominant in the transverse 
response. 



w 2 A xx 



co 2 A xz = oi 2 A M 



1 + 



fik 2 c 2 



■<U&(l+&r) 



n kc\ k^k. , , , 
, \2-i- — \-^- +k x k-c 2 
k 2 Al\ 4 to I k 2 



ci) - k z v b 
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pk 2 c : 



+ 1 4 yU kc 



A: 2 c 2 



10 — k z Vb 



(w - fe z u b ) 



An important observation is that <y p b < w p i, so that one can neglect 
\u)\ 2 in front ar pi . Furthermore, at equipartition, the ions contribute 
strongly to the instability and the typical wavenumber k x ~ to pl , 
whereas for p <k 1, the response of the electrons dominate, and 



A' 



the latter corresponding to the relativistic elec- 
tron plasma frequency. Thus the dispersion relation can be written 
(omitting also a term in f 2 r , and assuming klS 2 <K 1): 

/ , >2 



1 



pk 2 C 2 



4 /xkc 



k~c~ 2 k z i>b 
Sw 2 y 2 So) 



+2f k '~ Vh [k 2 S 2 i- - ^ 
So) \ 1 4 pkc k 2 
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We solve that equation by setting So) = v + iy with the approxi- 
mation |v| <K \y\ (y should not be confused with Lorentz factors 
appearing elsewhere). The growth rate is the positive root of the 
equation 

3 ° 

n y 



4 p&c 
which gives 



3+(l+^) 



7 

k 2 c 2 



. 



£ a u b o)„ 



\l + (k x S,) 2 ] 112 



(A9) 



(A10) 



provided (k x Sj) 3 > yJ^Zlp. If p > V?cr, the latter inequality is 
verified for k x Sj > 1, which means the filamentation instability is 
of the ion-ion type, with typical wavenumber k x ~ Sj l . If, however, 
fi < V?ci-> the instability is pushed towards higher values of k x , as 
the response of the hot electrons dominates. 

The frequency shift is obtained at the first order of the expan- 
sion in v: 

f y + (kSCPPSB** (AU) 



-k z v b 



3 fy + (k x S,)[l+(k x S,) 2 ]fio) pl 



Using the value of y and simplifying 1 + (k x 8f) 2 ~ (k x Sj) 2 , one 
obtains 



-k-v b -JJ^ 



(A12) 



which takes different scalings, depending on the comparison be- 
tween p and £ cr . 

If p > V?cn the filamentation instability grows at k x Sj ~ 1, so 

that 



-k z u b 



(p > VS. ~ i) 



If yt < V?ci- 5 maximum growth takes place at k x Sj - /j 



(A13) 
" 2 , so 



that 



v ~ -k.v b + fc) {fi < VS. k x S t ~ fi- [/2 ) , (A14) 

although it should be noted that the growth rate increases weakly 
with k x Sj, and that the frequency shift evolves in a non-trivial way 



with k x in the interval gj 6 p l/3 co p 



■ A' 



-I/ 2 /, 



We expect that the size of the Weibel filaments, i.e. the trans- 
verse coherence length is determined by the maximum k x ; thus 
£± ~ Sj far from equipartition, with t ± — » Sj close to equipar- 
tition. For this typical size, the filament Lorentz factor goes from 
7m ~ S" 2 when p < £ cr to y m ~ {p£ a y" 4 when % a < fi < 
and to y m ~ t; a l,4 p [/2 for p > ^|^. 

The wave-length in the normal direction is limited by k z v b < y 
which implies C\\ > SJ V^, which is comparable with the growth 
length. 



APPENDIX B: CO- VARIANT TRANSPORT EQUATION 
FOR PLASMA MODES 

We consider the transport equation of plasma modes with the right 
hand side involving the Einstein coefficients a and fi, knowing that 
in case of an instability, /3 = -2yj nst . The evolution is described 
through the occupation number per state, N(k, x), the energy den- 
sity e w being obtained by 



tia) k N(k, x) - 



£k 



and the energy flux S 



J %co k v g 



(k)N(k, x)- 



d 3 k 



(Bl) 



(B2) 



'{2nf 

The evolution equation is thus 

^-N(k,x) + v g • ^-N{k,x) = a{k) - j3(k)N{k,x) 
at ox 

we introduce the Lorentz factor y g associated with the group ve- 
locity and define the unitary 4-vector (u^) = (y g ,y g v g ). Then the 



(B3) 



(B4) 



equation can be rewritten as 

u^,N = y g (a-j3N)- 

The left hand side of the equation is clearly co- variant. Since N is 
invariant, the products y g a and yJ3 must be invariant too. In partic- 
ular we have the important property that the growth rate yi nst of an 
instability is transformed such that the product y g y ms t is invariant. 
This invariant product can be interpreted as the growth rate mea- 
sured in the wave-frame. 

Under a Lorentz transform, the energy density is transformed 
as follows: 

C , , , , d 3 k' r , , , , , d 3 k' 

e' = I ho)[N'(k',x') =■ = y'Jui)' t N'(V,x') T 

J k \2n) 3 J /s " y' g (27x) 3 

, / v g -v \l k-v \ , d 3 k 

\ c 2 /\ (o k j y,12nf 

I 



f(k)fUo k N(k, x) 



d 3 k 



(B5) 



where the factor f(k) = (l - ^)(l - ^). When the phase 



velocity is aligned with the group velocity and v,/,u g 



then 



f(k) = (1 - v g ■ v /c 2 ) 2 and y 2 J(k) = y g . 

And the energy flux is transformed as follows: 

For Weibel modes, that are not eigen-modes of the plasma, it is 

not clear that such description applies. However this is somewhat 

simpler, in the case we are treated in this paper, for they have a 

phase velocity and a group velocity equals to 



v m (k ± ) = v b \l -£ r /i 



2 + k\S 2 



l+k 2 S 2 



and thus the energy spectrum density 5 (k ± ,x) evolves according to 
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d 

-S + v m ■ VS = 2y inst S , (B6) 
at 

and multiplying both sides of the equation by y ra , the Lorentz factor 
associated with v m , we obtain a co-variant equation: 

t/^V„S = 2 rmri „ st S , (B7) 

where (/„, is the unitary 4-vector associated with the motion of 
the magnetic disturbances. Thus, without knowing how S is trans- 
formed by Lorentz transformations, we can conclude again that the 
product y m 7i nst is invariant. 



